Exact Tail Asymptotics of Dirichlet Distributions 



Enkelejd Hashorva 

Department of Mathematical Statistics and Actuarial Science 
University of Bern, Sidlcrstrassc 5,3012 Bern, Switzerland 

April 20, 2010 

, Abstract: Let X be a generalised symmetrised Dirichlet random vector inlR!',k > 2, and let tt„ GlR'',n > 1 be such 

, that lim„_).oo -P{X > Un} = 0. In this paper we derive an exact asymptotic expansion of P{X > u„} as n — > oo, 

, assuming that the associated random radius of X has distribution function in the Gumbel max-domain of attraction. 

, AMS 2000 subject classification: Primary 60F05; Secondary 60G70. 

, Key words and phrases: Dirichlet distributions; Elliptical distributions; Gumbel max-domain of attraction; Quadratic 

, programming; Exact tail asymptotics. 

; 1 Introduction 

, Let X = (^1, . • . , Xk)^ be a random vector \nF& , k > 2, and let m„, n > 1 be a positive sequence converging to infinity 

, as n — > oo. In our notation ^ stands for the transpose sign, and = below for equality of the distribution functions. For 

, any vector a = (ai, . . . , Ofc)^ G J?/^ \ (— oo, O]'^ the events 

! {X > u,ia} := {Xi > u^ai,. . .,Xk > u^ak}, n>l 

, are absorbing, i.e., lim„_j.oo P{X > M„a} = 0. For such instances it is of interest to determine the rate of convergence to 

[ of P{X > Una} as n — > oo. If X is a standard Gaussian random vector \nF& with non-singular covariance matrix E, 

] then we have the stochastic representation (see e.g. Cambanis et al. (1981)) 

; X^RA^U, (1.1) 

' with i? > such that R^ is Chi-squared distributed with k degrees of freedom, A a square matrix satisfying A^ A = S, 
' and U ~ {Ui, . . . , Uk)'^ a random vector uniformly distributed on the unit sphere oilR^ being independent of R. 

Results on the tail asymptotics of Gaussian random vectors are well-known, see e.g., Berman (1962), Dai and Mukherjea 

(2001), Hashorva (2005), Lu and Li (2009) and the references therein. 

Indeed, the radial decomposition in (|1.1[) is quite crucial also in an asymptotic context; it allows us to consider a general 
random variable R with some unspecified distribution function F . For such instances X is an elliptical random vector. 
If F is in the Gumbel max-domain of attraction (see (|3.13p below), then Theorem 3.1 in Hashorva (2007) implies 

P{X>Una] = {l + o{l))^{un)P{R> ^lUn}, n^oo (1.2) 

for any a ^IR^ \ (— oo, O]'^, with ^E* a known function and = (a^E^^a)^/^ with a the unique solution of the quadratic 
programming problem 

V{Tj^^,a) : minimise x^T,~^x under the constraint x £ [ai,oo) x • • • x [0^,00). (1-3) 

In this paper we are interested in extending (|1.2p for X with stochastic representation ()1.1|) where is a symmetrised 
Dirichlet random vector with parameter a G (0, 00)*"'. In the literature such X is referred to as a generalised symmetrised 
Dirichlet random vector, introduced and discussed in detail in Fang and Fang (1990). 
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Our novel contribution here is the derivation of the asymptotic expansion of P{X > UnO,}, 71 — > oo which can be described 
by (|1.2p also in the more general setup of Dirichlet random vectors. The main difference to the elliptical case (retrieved 
for a with components equal 1/2) is that some more general function (instead of appears in the asymptotics. 
Apparently, our new result cannot be derived by the existing one mentioned in (jl.2[) . since the new function ^E'q, depends 
explicitely on ex. Surprisingly, the constant fi does not depend on a. 

Another remarkable fact is that the Dirichlet random vector X and the associated elliptical random vector X* = ARV, 
with V uniformly distributed on the unit sphere oilR^ being independent of R can have the same tail asymptotic behaviour 
(up to some constant) if the index set {i : {Aa)i = 0, 7^ 1/2, 1 < i < fc} is empty. 

To this end wc mention that some possible results and application related (|1.2[) are: a) density approximation for Dirichlet 
random vectors, b) asymptotic expansions of conditional distributions, c) identification of the distribution of minima of 
Dirichlet random vectors, d) asymptotics of concomitants of order statistics, e) estimation of conditional distribution and 
conditional quantilc function, and f) asymptotic independence of sample maxima. In this paper we present few details 
regarding the last two applications in Example 2. 

Organisation of the paper: In Section 2 we give some preliminaries. The main result is presented in Section 3. We provide 
two illustrating examples in Section 4. Proofs are relegated to Section 5 followed by an Appendix. 



2 Preliminaries 



We shall introduce first some notation. Let in the following /, J be two non-empty disjoint index sets such that / U 



J = {1, . . . ,k},k > 2, and define for x ~ (xi,.. 



G -2?''' the subvector of x with respect to the index set / by 



xi := (xi, i E /)^ eF^ . If S eM''^'' is a square matrix, then the matrix S/j is obtained by retaining both the rows and 
the columns of S with indices in / and in J, respectively. Similarly we define Sj/, Sjj, S//. For notational simplicity we 
write xj instead of {xj)^ , (Sj,/)^-'^, respectively. Given x,y eM^ wc define 



X 
X 

x + y 
cx 
xy 

II^/IP 



> y, if Xi> yi, Vi = 1, . . . , fc, 

> y, if > T/j, Vi = 1,...,A;, 
= {xi+yi,...,Xk + ykV, 

= (cXi, . . . ,CXfc)^, ceR, 

= {xiyi, . . . ,XkykV , x/y := {xi/yi, . . . ,Xk/ykV ■ 

= xJY/jj Xj. 



Further, we write (and not e^./) for the ith unit vector in J?'^', where |/| > 1 denotes the number of elements of /. 
We shall be denoting by Ba.b a Beta random variable with positive parameters a and h and density function 



'{1-x) 



b-1 



xg (0,1), 



r(a)r(6) 

with r(-) the Gamma function. Further we write Y ^ H \i the random vector Y E]R^,k > 1 has distribution function 
H, and set H =1- H. 

Throughout this paper a := (ai, . . . ,afe)^ stands for a vector inM'^ with positive components, and 

k 



ax 



i^c{i,...,fc}. 



When K is empty, then ax equals 1. 



Definition 2.1. A random vector 14 = {Ui, . . . ,lAk)^ inW& , k > 2 is said to have symmetrised Dirichlet distribution with 
parameter a (henceforth 14 ^ S'D{k, a) ) if 14^14 = 1 almost surely, and {Ui, . . . ^lAk-i)^ possesses the density function h 
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given by 

M-i,. := 7^I^^(l-E«0""nV.p--\ E-Bl- (2-4) 
lli=i IQ^ij i=i i=i i=i 

Note that if a = 1/2 with 1 := (1,...,!)^ eM'', then S'D{k, a) is uniformly distributed on the unit sphere oiM^ . 

Definition 2.2. A random vector X inW&,k > 2 is said to possess a generalised symmetrised Dirichlet distribution 
(QSD) if X = ELI, where R > almost surely being independent ofU, and R ^ F,U ^ S'D{k,<y.). Write next 
X - gSV{k,a,F). 

In the fohowing we focus on distribution functions F with an infinite upper endpoint. Referring to Fang and Fang (1990) 
the density function (when it exists) of a QSV random vector can be defined via a density generator g and the parameter 
vector a e M'' . Spccificahy, let 5 be a positive measurable function such that 

nlir(a,) r _ 1 ^.-^a, (2.5) 



f g(x)x--idx = l, a:=5:< 

r(a) Jo ^ 



holds. In view of the aforementioned paper, see also Theorem 1 in Hashorva et al. (2007), if X ^ QS'D{k, a, F), then X 
possesses the density function h with density generator g defined by 

k k 

h{x) -.^ gC^\x,\^)Y[\x,\^"^-\ \fx,elR,l<t<k, (2.6) 
if and only if F possesses a density function / given in terms of the density generator g by 

m = 2nk£(^5(,2),2^-i^ Vre(0,oo). (2.7) 

r(a) 



Any subvector of a QSV random vector is again a QSV random vector (see Lemma [6731 in Appendix). 

A canonical example of a QST) random vector is a Kotz Type I QSD random vector in 5?'^ with density function 

{N-\-a)/s -n{ — \ ^ 

sT{{N + a)/s) Y[;^^ T{a,) fj[ 
with r > 0, s > 0, iV > —a, and density generator g given by 

g{x) = ccc^ exp(— ra;"), a;>0,c>0. 

In the standardised case = and 2r = s = 1 the random vector X possesses independent components with 

|A,p - GamTOa(a,, 1/2), Vi = l,...,fc, (2.9) 

where Gamma{ai, 1/2) is the Gamma distribution with parameters a^, 1/2. 

If A eJR^^^ is a non-singular square matrix, then A^ X possesses the density function 

hAix):=(iy ^^"^ exp(-||a:||V2)ni(Ca:).|^"'-\ WxeM'^, (2.10) 

V^/ nj=ir(a,;)|S|i/2 

with C := (A^^)^ , S := A'^ A. In particular, if ai = 1/2, 1 < i < k, then Ha is the density function of a centered Gaussian 
random vector with covariance matrix S. 



4 



3 Main Result 

Consider a QSV random vector X inM'^ , k > 2 with stochastic representation 

X^A^RU, (3.11) 

where R ^ F, k > 2 is independent oiU ^ S'D{k, a), a e (0, 00)*^', and A is a non-singular fc-dimensional square matrix. 
Without loss of generality we assume in the sequel that S := A is a correlation matrix, i.e., all the entries of the main 
diagonal of S are equal 1. 

For a given sequence u„a, n > 1 of thresholds inJ?'^ we are interested in the asymptotic behaviour of the joint survivor 
probability P{X > Una}. As in the elliptical setup (see Hashorva (2007)) it turns out that the tail asymptotics under 
consideration is closely related to the solution of the quadratic programming problem (jl.3p . If 



E-^a > 0, := (0,...,0)^ GiR'^ (3.12) 

is satisfied, then the minimum of the quadratic programming problem (|1.3p is attained at a, otherwise there exists a 
unique non-empty index set / C {l,...,fc} which defines the unique solution a of 7^(1]^^, a), see Proposition 16.21 in 
Appendix. In the following we refer to the index set / as the minimal index set. 

The only asymptotic assumption imposed below is that F is in the Gumbel max-domain of attraction with some scaling 
function w (for short F G GMDA{w)), i.e., 

F(u + x/w(u)) . ^ ,, „ , ^ 

lim ^ L ^ ^ " = cxp(-a;), VxGJ?, (3.13) 

u^oo F{u) 

with F -.= 1- F. See Reiss (1989), Embrechts et al. (1997), Falk et al. (2004), De Haan and Fcrreira (2006), or Resnick 
(2008) for more details on the max-domain of attractions. 

In order to avoid repetition we formulate the following assumption on the distribution function F. 

Assumption Al. is a univariate distribution with infinite upper endpoint such that F{Q) ~ satisfying p.l3p . 

In the sequel u„, n > 1 is a sequence of constants converging to infinity, and a e M'^ \ (— cx), 0]*^ is a given vector. If p.l3p 

holds, then we set 

Sn ■■= Un\\ai\\, Cri w((5„), A„ := (5„Cn, n>l, (3.14) 

with / the minimal index set of 7^(1]"^, a). Wc note in passing that in view of Proposition 16.21 ||a/j| G (0, 00). 

Below we show that also the parameter a ol X plays a crucial role in the tail asymptotics of interests since the following 

two index sets 

i := {1 < i < fc : ^ 1/2, (Ca)i 0}, Af {1, . . . , fc} \ L (3.15) 

appear explicitly in our asymptotic expansion. It is surprising that when the index set L is empty, then the random 
vector X and the associated elliptical random vector X* = ARV, with V ^ SV^k, 1/2) independent of R have the same 
tail asymptotics (up to some constant), i.e., lim„__i.oo P{X > u„a}/P{X* > M„a} = /i G (0, 00). We state next our main 
result. 

Theorem 3.1. Let X be a Dirichlet random vector inM'^ , k > 2 defined in p. lip , and set C := {A^)~^ , S := A'^ A. For 
a given vector a ^ F& \ (— oo,0]'', let I with m elements be the minimal index set corresponding to 7^(2^^, a) with the 
unique solution a > a. Suppose that Assumption Al holds, and define (5„, C„, An as in (|3.14p . Let w„, n > 1 6e a sequence 
of thresholds in IR^ such that 

lim Cn ("n - r = g/Gl?™, (3.16) 

n— ^C30 ^ 
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and if m < k 



lim(^ («„-u„a)^ = q,,G [-00,00)'=-'", J := {1, . . . , fc} \ /. (3.17) 



n— >oo 



a) If J := {1, . . . , k} \ I is non-empty and L <Z J with L, M defined in p.lSp . then (set u a/||a/|| J 

r(a) exp(-g7l]7/Mj) 
2ntir(a,)|S|V2n,^^e7l]7/«, 

^;^l-|/|-|J|/2+|L|/2-a.;p(^^)^ n ^ oo, (3.18) 



P{X > Un} = (1 + 0(l))r/ LC*M 1 r 



where 



Tj,L ■■= 



'M ■- 



I \{\{CjjyM^'''-'eM'y]{^~^)jjyjmdVj e (0,oo), 



b) If J is empty, then as n —>■ oo 

P{X>«4 = (1 +o(l))r,rX pr\|^n// »"''"'^'"'"^^('^")^ (3-19) 

:=4>,n.6Ll(Cy)«P"--^cxp(-MTs-iy)dye (0,oo). 



where 



In our notation JliGK ^« ^'^^ ^ — "^li^n K is empty. Next, given a G (0,00)*^ we define another vector 
a G (0,00)^= such that otM := 1a//2, and if L is non-empty set otL ckl. 

Corollary 3.2. Under the assumptions of Theorem \3.1\ (and the same notation), if further J is non-empty and either L 
is empty or 

Cj)iCjjV = (S-^)jj, (3.20) 

then we have 

P{X>u.} . (1 + o(l))rl,2-.--^ eM^,J^u'u,)P{Yj > g,} ^ £(1/2) 

xAi-l^l-l-^l/2+|i|/2-".F(,5„), n^oc, (3.21) 

where Yj = AjjY*j with Yj a standard Kotz Type I QST> random vector with parameter otj. Set P{Y j > q j} to 1 if 
|J|=0. 

Remarks: (a) By the properties of the scaling function w we have (see e.g., Resnick (2008)) 

lim uw{u) = 00. (3.22) 

U— f CO 

(h) If F is a univariate distribution function with upper endpoint 00, and further F G GAIDA{w), then for any r G 
(l,oo),r/ GiR 

^.^(^M^2)!ZM ^ Q ^323) 

x-i.00 F{x) 

holds, see Appendix Al. 

(c) Condition p.20p is satisfied in the special case J has only one element and T, is a correlation matrix. 

(d) In the 2-dimensional setup if A = (aij) j,j=i,2 is given by 

ail = 022 = 1, ai2 = a ^ 0, 021 = p ^ a, 
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then clearly A is non-singular and C :~ (A^) ^ has elements 

Cii = 022 = 1/(1 - crp), Ci2 = -/9/(l - crp), C21 = -cr/(l - CT/Cl). 

If a (l,a)^,a G (— oo,l], then the index set L in Theorem \3.1\ is non-empty for a = l/a,a ^ and ai ^ 1/2. With 
this choice of the constant a it follows that also the index set J is non-empty. Consequently, for the bivariate setup it is 
not possibly to have L non-empty and J empty. 



4 Examples 

We illustrate our result with two examples. First we consider the multivariate setup choosing the parameter vector a to 
have identical components, and then we deal with the bivariate setup. 

Example 4.1. Let X be a k-dimensional random vector as in Theorem \S.l[ where a ~ pi. ,p € (0, oo). We suppose 
that the matrix A is such that S = AJ A is given by 

E = pll' + {l-p)E, pe (-l/(fc-l),l), (4.24) 

with E ^W&^^ the identity matrix. 

Let Un u„l,n > 1 with Un,n > 1 given constants converging to infinity. In view of our asymptotical result we are 
able to derive the asymptotic of P{X > as n — > oo. We consider first the quadratic programming problem V{T., 1). 

Since the inverse matrix of S is 



l-p (l-p)(l + (fc-l)p) 
we obtain 



Bl . {pll- + il-p)Eiy\^ ^^ ' >0. 



{k-l)p 

Consequently, condition p.l2p holds implying that the unique solution o/7'(E,l) is 1 with the minimal index set I = 
{1, . . . , fc}. If p ^ 1/2, then the index set L is empty, and X is an elliptical random vector. We consider next the case 
p ^ 1/2 implying that L is defined by 

L:={i: {Cl)i = 0, l<i<k}, with C := (A^^)^. 

Suppose that C is a lower triangular matrix. Since B = C is positive definite, the matrix C can be explicitly determined 
by the well-known Cholesky decomposition. It follows easily that Ckk = Vbkk and 

Cfci = Cfc2 = • • • = Ck,k-i = bki/ \fbkk- 

Proceeding analogously we obtain 

Cii = = ■ ■ ■ = Cis-i, Cii + (i - l)c,i > 0, i = 1, . . . , fc - 1. 
Hence the index set L is empty, and 

\(Cl)i\ = cu + {k- i)cii, i = l,...,k. 

Since further 

llif = i^s-ii = ^^^^-^, /^=(i-p)('^-i)/vi + (fc-i)p, 
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p.l9p implies as n —> oo 



11 (^^^ r(M[fc/(i + (fc - i)^)]^-'^/^-^) U^=,\cu + ik- ^)c^l?''' . 



where 

k 



When p = 1/2 ('JC &em(; t/ius elliptically distributed) we obtain 



Vn:={un\\l\\w{un\\l\\Y-\ n>l. 



P{X>«„1} = (l + o(l)) ^ ^J »„u;KHlH))^-''fKHlH), ,woo. 



27r'^V2y|s| (1 ' S-iei)^ 

#o<e i/iai in the Gaussian case (i.e., is Chi-squared distributed with k-degrees of freedom) we have 



F{un) = (1 + o(l)) ^fc/2-rr(fc/2) exp(-</2), n ^ oo. 

Example 4.2. Consider X = RA^U a bivariate Dirichlet random vector with R ^ F and U ^ S'D{{ai, 02), 2). Assume 
that F satisfies the Assumption Al, and define the matrix A by 

an = 1, ai2 = p, 021 = 0, 022 = v^l - p"^, 

Consequently, S := A'^ A eM^^^ has elements 

Cll = (722 = 1, Cri2 = (721 = p 

and C := (A^)^^ /las elements 

cii = 1, C12 = 0, C21 = -p/\/T~f?, C22 = l/VT~~p2 



focus next on the asymptotics of P{Xi > u„,X2 > aii„},n — > 00 where a £ (— 00, 1] and w„,ri > 1 is a positive 
sequence of constants converging to infinity. Depending on the constant a we need to consider three cases: 

Case p < a: R follows easily that a = (l,a)^ is the solution of 'P{T,~^ ,a) with minimal index set I ~ {1,2}. Further, 
both J and L are empty and 

Sn = Uniil-2pa + a'^)/{l-p^)y/^ ^cu„, c ||a|| > 0, n > 1, 

_ C^{1-pY , ._ 2-2a,-2a, ( ^'P \ 

{l-ap){a-py \^T^) 
Consequently, in view of Theorem \3.1\ we obtain 

x{UnW{cUn))~^F{cUn), n ^ OO. 



In the special case a = 1 we have c = \/2/{l + p), and we may further write 

P{Xi>«„,X2>U„} = (1 + o(l)) J^^"^ + (1 - pr-\l + p)2-o.^3^2.,^2a. 

2r(ai)r(a2) 

x{UnW{cUn)Y^F{cUn), U^CO. 

Let bni '■= H~^{1 — l/n),n > 1 with Hy^,i = 1,2 the generalised inverse of the distribution function of Xi. R follows 
further that both Xi and X2 are in the Gumbel max- domain of attraction with the same scaling function w. In view of 
p.23p and the above asymptotics 

lim nP{Xi > bni,X2 > 6„2} = 0, (4.25) 
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hence Xi , X2 are asymptotically independent. 

We note in passing that the asymptotic independence is a crucial property of sample extremes with certain consequences 
for statistical applications. See Reiss and Thomas (2007), Hiisler and Li (2009), Das and Resnick (2009) and Peng 
(2010) for recent developments and applications. 

Case p = a Cz (—1, 1)- The minimal index set is I := {1} and if a2 7^ 0.5, then 

J = L = {2}, and||a/||=l, (Ca)i = 1, (Ca)2 = 0. 
In view of Theorem \3.1\ for any q2 € (— oo,c)o) we obtain thus (set r„ := Unw{un)) 

P{Xi> Un,X2> pUn + qUn/Vn} = (1 + o(l)) 



r(ai) 



where Y is symmetric about such that ^ Gamma{a2, 1/2). If 012 = 1/2, then L is empty. Also in this case the above 
asymptotics holds. Since Xi ^ R^Z, with Z = Bai^a2 being independent of R, in view of Theorem 12.3.1 of Berman 
(1992) we obtain the following conditional limit result 



lim P{X2 > pun + x^yun/w{un)\Xi > u„} = P{^/l~^Y > x} , ^fx^IR. (4.26) 

n—^oc 

Case a G (—00, 1), p £ (a, 1); The only difference to the above case is that for the choice of the threshold Un = Un(l, a)^ 



the vector q as defined in Theorem \3.1\ has components qi — and q2 ~ —00. Hence we have 

P{Xi>Un,X2>aun} = (l + o(l)) ^) \ ^ ' r''^^ F (Un) , n^oo. 

r(ai) 

5 Proofs 

Proof of Theorem 13.11 The idea of the proof is the same as that for the ehiptical setup (Theorem 3.1 in Hashorva 
(2007)), therefore we give the sketch of the proof omitting several details. 

Let Y* ^ QST){k + 1, a*, G) be a fc + 1 dimensional random vector where a.* — (ai, . . . , a^, 1)^ and distribution function 
G has the asymptotic behaviour 

G{u) = {l + o{l))uw{u)F{u) J^}^] u^oo. (5.27) 

2i (a + 1) 

By Lemma 7.4 in Hashorva (2007) such G exists. In view of (|6.40p the density function of I' = A^Y*^, N :— {1, . . . , /c} 
is given by 

where G := {A^)^^, S := A^ A, and ||a;|| := x^T^^^x, x £W& . By the properties of QST) random vectors and Lemma [6T] 
it follows that 

P{X>u,,] = (l + o(l))P{Y>M„}, n^oo. (5.28) 

Let / be the minimal index set of the quadratic programming problem P(I]^^, a) with the unique solution a. Applying 
Proposition 16 . 21 we obtain 

ai = ai, \\ai\\ > 0, aj = Y.jiT,Jjai > aj. 

Set for n > 1 

S„:^Un\\ai\\, C„ := w(u„j|a/i|), ^„ := w„||a/||ii;('u„||a/||) 
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and define v„ £ M!^ by 



(t>„)/ := t)„,/ = Cril/, {Vn)j -^Vnj ^[ — \ 1„ 



Cn 



5n 



1/2 



5n 



By Proposition 16 . 21 for any y £]R (set u := a/||a/||) 



|u„a + y/^;„|| = (5„ + (1 +o(l))[M7s^/yj + yJ(S ^)jjyj/2], n ^ oo. 



Further, p.22p and L d J with J non-empty imply 



ai,-|L|/2 



The above asymptotic might not hold in a subset ofF& with Lebesgue measure 0. If the index set J is empty and L is 
non-empty we have a different result, namely 

fc 



i=l ieM 

We consider first the case J is non-empty. For notational simplicity define 



ieM 



Next, for all n large we obtain 
P{Y > w„} 



r(s + i) 



dx 



r(a + l) 



ntir(a.)|s|V2 



r|,/|/2^-|/|-|,7|/2 



r(a + i) 



/-OO 

ni(CKa + yMO)«P"'"' / 



z-'" dG(z)/G((5„) 



Jy>v„{u„-u„a) -^j^ Jq{yj) 



dy, 



with q{yj) := (1 + o(l))[M7Sj/y7 + yJ(S ^)jjyj/2]. By the assumptions ujEj^yj + yJj{T, ^)jjyj/2 is positive for 
any y G-K*'' such that yj > 0/ and further the components of 'EJ^uj are all positive. The tail asymptotics of G and the 
fact that w is self- neglecting i.e., (see e.g., Resnick (2008)) 

w(u + x/w(u))/w{u) — > 1, li — !■ oo 



locally uniformly inK, imply G £ GAIDA{w), hence 

G{6n + t/U-G{6n + s/Cn) 



lim 

n— >oo 



exp(— s) — exp(— t), t > s, t, s S IR. 



Since further 



G(<5„) 

lim SnCn = oo, hm t»„(M„ - u„a) q 
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Fatou's Lemma and (|5.27|) yield 

. ^ P{X > n4 
limmi = 

T(a + 1) , 



X lim inf 

n— >oo 



/ _ r I[\(^Jjyjy^\^"'~'(^n + s/Cn)-^"dGiS„+s/Cn)/GiSn) 



dy 



r(a + 1)t* 77„(5-2" , T 1 s 



2nlir(aO|Sp/2n,,,e7E7>/"''"^ G(^„) ' 

with ■.= J^^^^^U,eL\iCjjyjh\''^^-'eM~yJji^~^)jjyj/2)dyj € (0,oo). 

The proof for the hmsup follows along the lines of the proof of Theorem 3.1 in Hashorva (2007), therefore omitted here. 
Hence using (|5.28p we have 

P{X>Un} = (l + o(l)) r ^ exp{-qjEjj^Ui)Bq, 

Similarly, if J is empty we obtain utilising further (|5.29p 

p{x>u^} = (1 + 0(1)) /^"^^ / Y[\{Cy).\'-'-'eM-y^^-'y)dy 

x(<5„Cn)'-'+l^l~'""n^„), n^oo, 

and thus the result follows. □ 
Proof of Corollary 13.21 If Y j is a standard Kotz Type I QSV random vector in JR'^"'" with parameters 5j, then 
p.lOp and (j3.20p imply that the random vector B^Y j,B := (C'Jj)^^ possesses the density function 

hB{x) := — ^exp(-a;J(E-i),,ja3j/2)ni(a7ja3j).|2"-i, ^xj eJR^-™. (5.30) 

Consequently, since 5,; = 1/2, i < k,i ^ L the constant tj.l in Theorem 13 . 1 1 can be re- written as 

rj,L = / \{\{Cj.jyM^^^''eM-yl,{^-^)jjyjmdyj 

= / Y[\{Cjjyj)^^'-'eM-yJji^-').jjyjmdyj = ^Sf^fe^J'i^./ > qA- 

Jyj>qj ,ej )JJ\ 
Hence we may further write 



P{X > u„} 
= (l + o(l))- 



r(a)i 



M 



2ntrr(a,)|S|V2n,,,e7S7/«, ' -^''^ KS-).,/^/^ 

x(5„c«)'-|'|-|''|/'-^^+|^|/m) 

^ n,e/r(«.)|S7/r/^n,e/-7S7/-/ r(a,) 

X (5„CJ'"l'l"l^l/'+l^l/'-^^ n'5„), n ^ oo, 
consequently the result follows. □ 
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6 Appendix 

Al. Proof of (|3.23p . A direct proof is given in Lemma 4.2 of Hashorva and Pakes (2009). Another (simpler) proof can be 
obtained by borrowing the argmnents of Lemma 2.2(b) in Resnick (2008). With Resnicks's idea it follows that for some 
e e (0, 1) small enough 

F{x + s/w{x)) ^ c 

F{x) - (l+esy-' ^ ■ ^ 

is valid for all s > and any x large with c some positive constant. p.23p follows now easily. 

A2. The next lemma is a simple generalisation of Lemma 7.3 in Hashorva (2007). 
Lemma 6.1. Let R, R* be two independent positive random variables such that 

1™ ^"l = ce(0,oo). (6.32) 

u^oo P{R* >u} ^ ' ^ ^ 

Suppose that R has the distribution function F satisfying the Assumption Al, and let 14 ^ S'D{k,cy.) being independent 
of R and R* . Then we have 

P{RAU > ua} , 
u^oa P{R*AU > ua\ 

for any a \ (— oo, 0]*^ and A £IR^^^ a non-singular matrix. 

Proof of Lemma [6. II Since R has distribution function in the Gumbcl max-domain of attraction, then by (|3.23p for any 
ii' > 1 we obtain 

P{R> Ku) , P{R* > Ku} 
lim ——7 ^ = Imi ——7 ^ = 0. 

u^oo P{R > u] u^x P{R* > u} 

This fact together with A being non-singular implies the proof. □ 
A3. The next proposition can be found in Hashorva (2005). 

Proposition 6.2. Let S G > 2, be a positive definite correlation matrix and let a E \ (— 00,0]*^ be a given 

vector. Then the quadratic programming problem 

V^S^^ ,a) : minimise \\x\\'^ under the linear constraint x> a 

has a unique solution a defined by a unique non-empty index set L C {1, . . . ,k} such that 

ai = ai > O7, Ef/a/ > O7, (6.34) 
min|la;|p = mina;^S"^a; = ||a/|p = ||a/|p = ajj:jjai > (6.35) 



and if \I\ < k, then with J :— {1, . . . , d} \ / 

aj - -((S-i)jj)-i(E-i)j/a7 = Sj/E7/aj>aj. (6.36) 

k 



Furthermore, for any x G IR we have 



aj^E^^a = x]T.-}ai = xjj:j^ai (6.37) 



and 2 < |/| < k, provided that a ~ cl, c G (0, oo). 



A4. The next result follows easily, see Fang and Fang (1990). 
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Lemma 6.3. Let X ^ QS'D{k, ot, F) be a random vector in]R^,k > 2 where F is a distribution function satisfying 
F{0) ~ 0. For any index I C {1, . . . , fc} with m elements we have 

Xi = RiUi, (6.38) 

with Ri being independent oflAi ^ S'D{k,cy.i), and Rj possesses the density function fi defined by 

fiiz) = 2z^"'- \ j _ r(r^ - ^2)«-o.-i^-2(o-i) ^ze{Q,^). (6.39) 

Furthermore for any non-singular matrix A G-ZR™^™ the random vector Xj possesses the density function Haj given 
by (set C :~ (A'^)"^ ,Yi := A^ A and write jSj for the determinant ofE) 

hAji^i) = ^'^^"^P'^^.'ff'y";"' . r (^^-|l^/rr"-V-^("-^)df(r), yxeRK (6.40) 

References 

[1] Bcrman, M.S. (1962) A law of large numbers for the maxiniuin in a stationary Gaussian sequence. Ann. Math. Stats. 
33,1, 93-97. 

[2] Berman, M.S. (1992) Sojourns and Extremes of Stochastic Processes, Wadsworth & Brooks/ Cole, Boston. 

[3] Cambanis, S., Huang, S., and Simons, G. (1981) On the theory of elliptically contoured distributions. J. Multivariate 
Anal. 11,3, 368-385. 

[4] Dai, M., and Mukherjea, A. (2001) Identification of the parameters of a multivariate normal vector by the distribution 
of the minimum. J. Theoretical Prob. 14,1, 267-298. 

[5] Das, B., and Resnick, S.I. (2009) Detecting a conditional extreme value model. Preprint. 

[6] De Haan, L., and Ferreira, A. (2006) Extreme Value Theory. An Introduction. Springer. 

[7] Embrechts, P., Klppelberg, C., and Mikosch, T. (1997) Modeling extremal events for finance and insurance. Berlin, 
Springer. 

[8] Fang, K.-T., and Fang, Bi-Qi. (1990) Generalised symmetrised Dirichlct distributions. In Statistical inference in 
elliptically contoured and related distributions, K.T. Fang and T.W. Anderson, eds, AUerton Press, New York, pp. 
127-136. 

[9] Falk M., Hiisler, J., and Reiss R.-D. (2004) Laws of Small Numbers: Extremes and Rare Events. DMV Seminar 23, 
Second edition, Birkhauser, Basel. 

[10] Fang, K.-T., Kotz, S., and Ng, K.-W. (1990) Symmetric Multivariate and Related Distributions. Chapman & Hall, 
London. 

[11] Hashorva, E. (2005) Asymptotics and bounds for multivariate Gaussian tails. J. Theoretical Prob. 18,1, 79-97. 

[12] Hashorva, E. (2007) Asymptotics properties of type I elliptical random vectors. Extremes, 10,4, 175-206. 

[13] Hashorva, E., Kotz, S., and Kumc, A. (2007) Xp-norm generalised symmetrised Dirichlct distributions. Albanian J. 
Math. 1, 31-56. 

[14] Hashorva, E., and Pakes A.G. (2009) Tail asymptotics under beta random scaling. Submitted to Extremes. 



13 



[15] Hiisler, J., and Li, D. (2009) Testing asymptotic independece in bivarite extremes. J. Stat. Plan. Inf. 139, 990-998. 

[16] Lu, D., and Li, W.V. (2009) A note on multivariate Gaussian estimates. J. Math. Anal. AppL, 354,2, 704-707. 

[17] Peng, L. (2010) A practical way for estimating tail dependence functions. Statistica Sinica, 20, 365-378. 

[18] Reiss, R-D. (1989) Approximate Distributions of Order Statistics: With Applications to Nonparametric Statistics. 
Springer, New York. 

[19] Reiss, R-D., and Thomas, M. (2007) Statistical Analysis of Extreme Values. From Insurance, Finance, Hydrology 
and Other Fields. Third Edition, Birkhauser, Basel. 

[20] Resnick, S.I. (2008) Extreme Values, Regular Variation and Point Processes. Soft cover edition. Springer, New York. 



